Estimation of parameters of gravitational waves from coalescing binaries 
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In this paper we deal with the measurement of the parameters of the gravitational wave signal 
emitted by a coalescing binary signal. We present the results of Monte Carlo simulations carried 
out for the case of the initial LIGO, incorporating the first post-Newtonian corrections into the 
waveform. Using the parameters so determined, we estimate the direction to the source. We stress 
the use of the time-of-coalescence rather than the time-of-arrival of the signal to determine the 
direction of the source. We show that this can considerably reduce the errors in the determination 
of the direction of the source. 



A major source for the planned interferometric gravitational wave detectors like LIGO ]l| and VIRGO |Q| is the 
radiation emitted by a binary system of stars during the last few minutes of its evolution just before the two stars 
coalesce. The interferometers are expected to have high enough sensitivity to detect sources at cosmological distances. 
Since the waveform emitted during a binary coalescence can be modeled fairly accurately it is possible to enhance the 
visibility or the 'signal-to-noise' ratic|^, (defined later in the text in eq. [l^) of the signal in noisy data by employing 
0^ , such powerful techniques such as Weiner filtering . With the aid of such methods one essentially compares signal 
energy to that of the noise as opposed to the case of a burst signal where one has to contend with the signal power 
vis-a-vis noise power [|| . The fact that the signal waveform can be predicted very well also implies that the parameters 
of the waveform can be estimated to a high accuracy. For instance, by tracking the radiation emitted by a neutron 
star-neutron starn (NS-NS) binary starting from 10 Hz. to say IkHz., the masses can be determined to better than a 
few percent [p|-[10[. Observation of several coalescence events could be potentially used to determine Hubble parameter 
00 to an accuracy of 10% Once in a while, when an event produces a very high SNR (> 40) it would be possible 
' to infer the presence of gravitational wave tails in the signal waveform and test general relativity in the strongly 
vn non-linear regime These are, but a modest list of interesting physical and astrophysical information that the 

observation of gravitational waves is expected to bring forth. 
O ■ In the recent past the problem of optimal detection and estimation of parameters of the binary has gained con- 
O^' siderable attention, thanks to the funding of the American LIGO and Franco-Italian VIRGO projects. The chirp 
^ [ waveform is expected to be detected with the aid of Weiner filtering, wherein the detector output is correlated with 
a set of templates which together span the relevant range of the parameters of the binary waveform. When a signal 
^ [ is present in the detector output, the filter matched to that particular signal will on the average obtain the largest 
' ^ ■ possible correlation with the detector output amongst all other filters: thus enabling detection as well as estimation 
' of parameters. The parameters of the template that obtains the maximum correlation is an un biased estimate of the 
actual signal parameters. Of course, the detector output also contains noise which can effect the correlation thereby 
giving rise to spurious maxima even when the parameters of the template and those of the signal are mismatched. 
Thus, the measured parameters are in error. These errors tend to reduce with increasing SNR. 

Based on analytical computation of the covariance matrix of the errors in the estimation of parameters and the 
covariances among them (for the advanced LIGO), it is now generally agreed that the chirp mass (which is a combi- 
nation of the masses of the binary stars) can be determined at a relative accuracy ~ at a SNR of 10. The 
reduced mass, as it turns out, will have a much larger error especially if the spins are large. In the computation of 
the covariance matrix one makes the crucial assumption of high SNR. As a result the covariance matrix only gives 
a bound on the errors in the estimation of parameters and a more rigorous or a different method of computing the 
variances and covariances is in order. Alternatively one can determine the errors through the Monte Carlo simulation 
method. The basic idea here is to mimic the actual detection problem on a computer by adding numerous realizations 
of the noise to the signal and then passing the resultant time series through a bank of Weiner filters each matched to a 
gravitational wave signal from a coalescing binary system with a particular set of parameters. The covariance matrix 



^henceforth denoted as SNR 

■^In this paper a neutron star is considered to have a mass of 1.4 M0,and a black hole will be considered to have a mass of 
10. Mq 
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can be computed by using the ensemble of measured values of the signal parameters. Such a method is robust in 
the sense that it does not make any assumptions about the SNR and therefore is much closer to the actual detection 
problem. It also squarely addresses the difficulties faced in a realistic data analysis exercise such as the finite sampling 
rate, the finiteness tof spacing in between the filters, etc. For details on the Monte Carlo method see Balasubramanian, 
Sathyaprakash and Dhurandhar 

In this paper we summarize the outcome, of the first in a series, of Monte Carlo simulations which suggests that 
the covariance matrix grossly underestimates the errors in the estimation of the parmeters by over a factor of 3 at an 
SNR of 10. 

The standard method to deduce the direction to the direction to the source of the gravitational radiation is to use 
the differences in the times-of-arrival in a network of three or more detectors. We point out that since the instant 
of coalescence of a binary system can be measured to a much greater accuracy than that of the time-of-arrival, the 
former can be employed to infer the direction of the source at a much lower uncertainty than has been so far thought. 

In the point mass approximation the restricted first-order post-Newtonian gravitational radiation emitted by a 
binary system of stars induces a strain in the detector which can be written as 
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where f(t) is the instantaneous gravitational wave frequency, which at this level of approximation is equal to twice 
the orbital frequency and is implicitly given by: 
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and ip(t) is the phase of the signal given by 
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where ta and <E> are, respectively, the so called time-of-arrival and phase at which the signal reaches a frequency fa] 
$ is a constant depending on the relative orientations of the binary orbit and the arms of the interferometer; tq and 
Ti are constants, having dimensions of time, depending on the masses mi and m2 of the binary system. They are 
referred to as chirp times and are given by. 
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where /i is the reduced mass, rj is the reduced mass divided by the total mass Af, and M = ^"^/^Np/^ is the chirp 
mass. In this level of approximation the post-Newtonian waveform is thus characterised by a set of five parameters: 
The amplitude parameter A, the time-of-arrival ta^ the phase at the time of arrival $, and the two masses mi and 
m2 . We shall find below a more convenient parameterization of the waveform. 

The Fourier transform of the waveform in the stationary phase approximation is given by 
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where A '\s a, normalization constant and is fixed by specifying the SNR, and tp/iif) are functions only of frequency 
given by. 
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Here are related to the set of parameters introduced earlier and are given by, 



A° =^ and A^ = {^c, *c, M, tj}, fi = 1, . . . ,4 , (U) 

where, 

tc=ta+To+ Ti (12) 

$c = $ - 27rfJa + ^^^y^ + (13) 

For / < 0, the Fourier transform can be obtained by the relation h{~f) = h*{f), obeyed by all real functions h{t). 

Central to the computation of the covariance matrix is the scalar product defined over the space of the signal 
waveforms. Given waveforms h(t] A) and g{t; A) their scalar product is defined by, 

oo _ 
fa 

where Sn{f) is the two sided detector noise power spectral density. The scalar product defines a norm on the vector 
space. The SNR can now be defined for a signal h{t; A) when it is passed through the matched filter as the norm of 
the signal i.e., 

p={Kh)^/\ (15) 
The covariance matrix C^^ is the inverse of the Fisher information matrix F^j, given by, 

r.. = (^,^); c,. = r-V- (i6) 

The diagonal elements of the covariance matrix provide us an estimate of the variances to be expected in the measured 
values of the parameters in a given experiment. In Fig. ^ we have shown the behaviour of the la uncertainties (square 
root of the variances) in the estimation of parameters ta,TQ,Ti and re, as a function of the SNR. As is well known the 
covariance matrix predicts that the errors fall off in inverse propotion to the SNR. It is to be noted thet this is valid 
only in the 'high' SNR limit. 

In the actual detection problem the detector output x{t) is filtered through a host of search templates corresponding 
to test parameters, (A and the template that obtains the maximum correlation with the output gives us the measured 
values mA of the signal. These measured values will in general differ from the actual signal parameters. With the aid of 
a large number of detectors one obtains an ensemble of measured values ^ A. The average of such an ensemble provides 
provides us with an estimate eA and the variances and covariances ax computed using ^A helps us in deducing the 
errors in the estimation and how the different parameters are correlated: 



eA^^A, ai = („,A-„,A)', S^'- ^ , (17) 

where an overbar denotes an average over an ensemble and the correlation coefficients. In reality we will have 
only a few detectors and hence a numerical simulation needs to be carried out to deduce the errors in the estimation 
of parameters. We have carried out such a numerical simulation by using in excess of 5000 realizations of detector 
noise. Thus, the results of our simulations are statistically significant, the errors in the estimation of various statistical 
quantities such as the mean and the variance are negligible and hence we do not plot the error bars in our curves. In 
Figure ^ the dotted line corresponds to the behaviour of the errors in the estimation of various parameters computed 
at several values of the SNR. At low values of the SNR (p ^ 10) there is a significant departure of the observed errors 
from those predicted by the covariance matrix. However at an SNR of > 30 the two curves merge together indicating 
the validity of the covariance matrix at high enough SNRs. Since most of the events which the interferometric 
detectors will observe are expected to have an SNR less than 10, we conclude that the accuracy in the determination 
of the parameters is not as high as was thought to be, based on the values of the covariance matrix. Detailed analysis 
suggests that this discrepancy is larger when higher post-Newtonian corrections are taken into account. Consequently 
a more exhaustive analysis than has been reported here or elsewhere is in order. We are in the process of carrying 
out simulations by taking a second order post-Newtonian waveform and including other physical effects such as the 
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With reference to Figure [T] we point out that the instant of coalescence tc can be determined to an accuracy much 
better than the time-of-arrival ta- TypicaUy ate ^ factor of 50 less than Ot^- Consequently with the aid of tc we 
can fix the direction to the source at a much higher accuracy than with ta as can be seen from equation (|T^), tc 
is the sum of tq, ti and ta, and the errors in these parameters tend to cancel because of the presence of negative 
covariances. It is to be noted, however, that as of now we do not know the orbit of the binary accurately enough, 
to predict the exact instant of the coalescence. In fact the frequency cutoff imposed by the onset of the plunge orbit 
will make it difhcult to calculate tc- Moreover, for a realistic detector the noise increases with the frequency beyond 
about 150 Hz. and the frequency which will be of more interest to us is the one where the power spectrum of the 
signal divided by the power spectrum of the noise in the detector reaches a maximum. This of course assumes that 
all the detectors to used for direction measurement are identical. 

Let t denote a convenient time parameter which if measured at three detectors gives the direction to the source. 
For our purpose we take this parameter to be either tc or ta- We will assume a very simple configuration of three 
identical detectors placed at the vertices of an equilateral triangle on the equitorial plane of the earth. We will take 
the separation between any two detectors as L and consequently the maximum time delay induced by i as A = i/c 
where c is the velocity of light. Figure |^ illustrates such a network and the coordinate system used. The X axis passes 
through the detectors labelled 1 and 2 and has its origin at the former. The positive Y axis is chosen as shown and 
the Z axis is perpendicular to this plane and forms a right handed coordinate system with the other axes. Given the 
values ti, t2 and as the measured values of t in the three detectors respectively we can deduce the direction to the 
source via the time delays "t — ^2 ~ ti and S — t^ — ti- As the noise in each detector is uncorrelated with the noise in 
the others ti , t2 and t^ are uncorrelated and is we assume that all the detectors are identical then the errors in t in 
all the detectors will be the same. However, 7 and S will have non-zero covariances. Thus, 
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The angles (f> and 9 are related to 7 and S by. 
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and, 
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The errors in the measurement of the time parameter will induce errors in the measurement of the angles. Assuming 
the errors in the time parameter to be small we can write the expressions for cta and ag as, 
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where, 
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The factors 5</,(0, 0) and ge{9, (f) are typically of order unity for most directions. For earth based detectors the value 
of A is « 15 ms. This certainly rules out the use of ta to determine the direction as the error in this parameter even 
for high SNRs is much more than A. Though for the initial LIGO we can use tc, the errors at an SNR of 10 will 
be around two degrees which is too large to make an optical identification of the source. In order to determine the 
direction to arcsecond resolution one needs the sensitivity of the advanced LIGO. 

We would like to end this paper with the remark that though the inclusion of higher order post-Newtonian corrections 
is expected to bring down the precision with which the astrophysical parameters can be determined, they are not a 
cause for worry for the detection problem. Even at the second post-Newtonian correction the signal can essentially be 
detected with the aid of a one-dimensional lattice of templates However, more extensive numerical Monte Carlo 
simulations need to be performed to gain further insight into the estimation problem. 
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FIG. 1. Dependence of the errors in the estimation of parameters of the post-Newtonian waveform i.e. {ctq , ''"ti , crt„ , crtp } 
as a function of SNR. The solid line represents the theoretically computed errors whereas the dotted line represents the errors 
obtained through Monte Carlo simulations. The simulations have been carried out for case of a IOM0 — I.4M0 binary system. 
The errors in the parameters are expressed in ms. 



FIG. 2. This figure illustrates the network of three detectors and the choice of the coordinate system employed. The Z axis 

is perpendicular to the plane of the paper and points upward. The angle 6 is defined as the angle which the direction vector 
makes with the positive Z axis and (p is defined on the X — Y plane as shown in the figure. 
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